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Abstract 

We observed the narrow-line quasar SDSS J094857. 3+002225, which has the highest known radio loud- 
ness for a narrow-line Seyfert 1 galaxy (NLS1), at 1.7-15.4 GHz with the Very Long Baseline Array (VLBA). 
This is the first very-long-baseline interfcrometry (VLBI) investigation for a radio-loud NLS1. We indepen- 
dently found very high brightness temperatures from (1) its compactness in a VLBA image and (2) flux 
variation among the VLBA observation, our other observations with the VLBA, and the Very Large 
Array (VLA). A Doppler factor larger than 2.7-5.5 was required to meet an intrinsic limit of brightness 
temperature in the rest frame. This is evidence for highly relativistic nonthermal jets in an NLS1. We 
suggest that the Doppler factor is one of the most crucial parameters determining the radio loudness of 
NLSls. The accretion disk of SDSS J094857. 3+002225 is probably in the very high state, rather than the 
high/soft state, by analogy with X-ray binaries with strong radio outbursts and supcrluminal jets such as 
GRS 1915+105. 
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1. Introduction 

Narrow-line Seyfert 1 galaxies (NLSls; Osterbrock & 
Poggc 1985; Poggc 2000) arc usually found as radio- 
quiet objects (Zhou & Wang 2002). Various properties 
observed through optical and X-ray bands suggest that 
NLSls are active galactic nuclei (AGNs) of highly ac- 
creting (Lboi/^Edd ~ 1' wnere iboi and Z/Edd are bolo- 
metric and Eddington luminosities, respectively) sys- 
tems (e.g., Pounds et al. 1995; Boroson 2002). The radio- 
quietness of NLSls may be related to the suppression of 
radio jets emanated from accretion disks with high ac- 
cretion rates (Greene et al. 2006), as well as the disks of 
X-ray binaries in the high/soft state (see, e.g., McClintock 
& Remillard 2003 for a review). 

On the other hand, radio-loud NLSls do ex- 
ist, although they are very rare (Zhou & Wang 
2002). There is no significant difference between 
radio-quiet and radio-loud NLSls in optical and 
X-ray properties (SDSS J094857.3+002225, Zhou 
et al. 2003; RX J0134.2-4258, Grupe et 



al. 2000; RGB J0044+193, Siebert et al. 
1999; PKS 2004-447, Oshlack et al. 2001; 
PKS 0558-504, Wang et al. 2001). The origin of 
the radio excess of the radio-loud NLSls has been 
unknown. Hardening of X-ray spectra during rapid X-ray 
flares of the radio-loud NLS1 PKS 0558-504 in a few 
minutes could arise from transient spectral dominance 
of synchrotron emission from relativistically boosted 
jets (Wang et al. 2001) similar to radio-loud quasars 
(e.g., Reeves et al. 1997). However, properties of radio 
jets in NLSls are slightly known. Their linear sizes are 
< 300 pc (Ulvestad et al. 1995), which correspond to 
~ 1 arcsec resolutions of connected interferometries; there 
are only a few images that have been obtained by very- 
long-baseline interferometry (VLBI) observations (Lai et 
al. 2004; Middelberg et al. 2004). 

We started VLBI-imaging several NLSls. 
SDSS J094857.3+002225 was discovered at z = 0.584 in 
the Sloan Digital Sky Survey (SDSS) early data release 
(Williams et al. 2002). With FWHM(H/3)«1500 km s~ 4 
and undetected [OIII] lines, this source was consistent 
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Table 1. Parameters of observations. 



Date 


Array 


v c 
(GHz) 


Au 
(MHz) 


pol. 


^scan X ^V scan 
(8) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


2003Oct30 


VLBA 


15.365 


32 


L 


46x52 






8.421 


32 


R 


35x18 






4.987 


32 


L 


56x7 






2.271 


32 


R 


85x6 






1.667 


32 


L 


94x6 


2005May09 


VLBA 


1.667 


32 


L 


46x12 


2003Decl5 


VLA-B 


14.94 


100 


dual 


43x2 






8.460 


100 


dual 


33x1 






4.860 


100 


dual 


33x1 






1.425 


100 


dual 


43x1 



Col. (1) observation date; Col. (2) array; Col. (3) center frequency; Col. (4) bandwidth; Col. (5) polarization. L and R denote left and 
right circular polarization, respectively; Col. (6) scan lengths X number of scans. 



with the conventional definition of NLSls (Zhou et al. 
2003). Its radio counterpart was found in the Faint 
Images of the Radio Sky at Twenty-centimeter (FIRST) 
survey (Becker et al. 1995). This has been the most 
radio-loud (R « 2000) object in NLSls (Zhou et al. 
2003), where R is "radio loudness" conventionally defined 
as the ratio of 6 cm radio to 5-band flux densities 
with a threshold of R = 10 separating radio-loud and 
radio-quiet objects (e.g., Visnovsky et al. 1992; Stocke et 
al. 1992; Kellermann et al. 1994). Zhou et al. (2003) sug- 
gested that its rapid (~years) radio-flux variations were 
caused by relativistic boosting of a jet. This is appropri- 
ate for the first target of VLBI observations, because of 
the brightest (~ 300 mJy) NLS1 in radio bands. In the 
present paper, we report VLBI observations in milliarc- 
second (mas) resolutions for SDSS J094857.3+002225; 
this is the first VLBI study for a radio- loud NLS1. 
Throughout this paper, a flat cosmology is assumed, 
with H = 71 km s" 1 Mpc" 1 , Q M = 0.27, and Q A = 0.73 
(Spergel et al. 2003). Therefore, 1 mas corresponds to 
6.6 pc at the distance to SDSS J094857.3+002225. Its 
optical magnitude ms(SDSS) w 18.89 mag (Zhou et al. 
2003) corresponding to Mb = —23.8 mag implies that this 
source is a quasar rather than a Seyfert galaxy according 
to the standard definition of a quasar, Mb < —23 mag. 

2. Observations and data reduction 

We observed SDSS J094857.3+002225 using ten sta- 
tions of the Very Long Baseline Array (VLBA) on 
October 30, 2003 and May 09, 2005. The observation pa- 
rameters arc shown in Table 1. Standard procedures were 
applied for data reduction using the Astronomical Image 
Processing System (AIPS; Greisen 2003). Amplitude cal- 
ibration used a priori gain values, together with system 
temperatures measured during the observations, and was 
typically 5% accurate. Imaging, deconvolution, and self- 
calibration for phase and amplitude were performed us- 
ing the Difmap software (Shepherd 1997). The final self- 
calibration in phase was done at all the frequencies with a 



solution interval of 10 s, which is much less than coherence 
times. 

SDSS J094857.3+002225 was also observed in a snap- 
shot mode with the Very Large Array (VLA) B-array con- 
figuration on December 15, 2003. The VLA observations 
were carried out in good weather conditions. Data were 
calibrated using the AIPS, according to the standard man- 
ner. Maps in Stokes I were made, and deconvolution and 
self-calibration were performed in the AIPS. A solution 
interval of 10 s, similar to the VLBA data reduction, was 
used in self-calibration to correct of visibility phases. The 
flux-scaling factors of visibility amplitude were determined 
by observations from the standard calibrator 3C 48, with 
an uncertainty of ~3%. 

3. RESULTS 

3.1. VLBA images 

In the VLBA data analyses, source structure models 
were established by visibility-based model fitting with 
point and/or elliptical-Gaussian profilc(s) in the Difmap. 
Here, we describe the results of the fitting, particularly at 
15 GHz, where the highest spatial resolution was available 
during the observations. The free parameters were flux 
density and two-dimensional position for a point source 
model, along with major axis, axis ratio, and position an- 
gle of the axis for an elliptical-Gaussian model. 

One point source resulted in x 2 /do,f = 77821/57109; 
adding one more point source showed x 2 /do/ = 
73870/57106: a dramatic improvement of A^ 2 /Adof = 
3950/3, which is significant at >99.9% confidence. We 
found both a strong component "CO" and a weak compo- 
nent "CI," not only at 15 GHz but also at 5 and 8.4 GHz. 
However, the two components were not separable from 
each other at 1.7 and 2.3 GHz, where the resolution was 
poorer (Table 2 and Fig. 1). No other significant emis- 
sion in the residual images was found. We had initially 
suspected that the weaker component was created by a 
deconvolution error, because the separation was the same 
as or smaller than the synthesized beams. However, a 
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Fig. 1. VLBA contour maps of SDSS J094857.3+002225 on Oct 30, 2003. The contour levels arc 3cr x -1, 1, 2, 4, 8, 16, 32, 64, 
128, and 256, where a is the RMS of image noise at each frequency. In the images at 5.0, 8.4, and 15.4 GHz, "+" symbols represent 
the positions of CO and CI components, which were identified by visibility-based model fitting (Section 3.1). 



number of test images from fake visibility data created by 
Monte-Carlo simulations with an identical (it, i>)-coverage 
and a similar sensitivity indicate that a strong compo- 
nent and a weak one that are separated by >0.10 mas at 
PA =30° can be clearly resolved from each other. We are 
confident, therefore, in the existence of the two compo- 
nents separated by 0.41 mas at 15 GHz. Although the 
separations between CO and CI were different at observ- 
ing frequencies, the position angles of CI from CO were 
consistent (Table 2). 

Furthermore, we attempted to fit the visibility data 
with two elliptical Gaussians. However, this attempt 
failed because the source sizes were too small. To evaluate 
how a large source can be resolved in the VLBA observa- 
tion, we made imaging simulations in the same manner 
using fake visibility data that were made by combining a 
strong circular-Gaussian source of 461 mJy with various 
diameters and a weak point source of 21 mJy separated 
by 0.4 mas at PA ~30° (see similar analyses in Taylor et 
al. 2004). These tests suggest that the stronger source 
can be clearly resolved if the given size is >0.05 mas in 
diameter. The VLBA observation should have the ability 
to resolve a size of ~0.10 mas for the stronger component. 
The size of 0.10 mas corresponds to a decrease of ~ 10% in 
visibility amplitude at the longest baseline of the VLBA 



at 15 GHz. Because of the sufficiently good sensitivity 
and calibration, we were able to determine the fringe vis- 
ibility to 5% (see similar discussion by Kellermann et al. 
1998). Thus, we assume the effective resolution of our 
VLBA observation at 15 GHz is around 0.10 mas, which 
is about quarter of the synthesized beam. The compo- 
nent CO is presumably compact with a size of < 0.10 mas, 
corresponding to <0.66 pc on a linear scale. 

3.2. VLA images 

In the VLA images, only a single, unresolved component 
was found at all the frequencies, except for 1.4 GHz, where 
radio interference affected the data. No diffuse emission 
component was found, even in (w,u)-tapcrcd images. Flux 
density was measured by Gaussian model fitting on the 
image domain using the task JMFIT of the AIPS. 

3. 3. Measured flux densities 

Measured flux densities are listed in Table 2. We 
estimated the total error from the root sum square of 
the uncertainty of amplitude calibration and the error in 
Gaussian fitting (including thermal noise). 

The total spectra of observed VLBA flux densities were 
slightly inverted (a > 0, S„ oc v +a ), which is consistent 
with past observations. Non-simultaneous observations 
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Table 2. Observational results. 



V 


Comp. 




s v 




<T rms 




^min 


PA 


Z(C0-C1) 


iM(CO-Cl) 


(GHz) 






(mJy) 


(mJy beam -1 ) 


(arcsec) 


(arcsec) 


(deg) 


(mas) 


(deg) 


(1) 


(2) 




(3) 




(4) 


* (5) ' 


* (6) ' 


(7) 


(8) 


(9) 


VLBA observed 


on Oct 30, 


2003 














15.365L 


CO 


461.3 


± 


23.1 


0.65 


0.0011 


0.0004 


-3.3 








CI 


20.6 


± 


1.6 










0.41±0.02 


+30.7 


8.421R 


CO 


363.0 


± 


18.2 


0.37 


0.0019 


0.0008 


-2.7 








CI 


25.2 


± 


1.4 










0.51±0.01 


+34.9 


4.987L 


CO 


258.9 


± 


13.0 


0.56 


0.0031 


0.0012 


-3.8 








CI 


32.0 


± 


2.4 










0.75±0.02 


+33.1 


2.271R 




245.4 


± 


12.3 


0.67 


0.0072 


0.0029 


-2.8 




. . . 


1.667L 




173.0 


± 


8.7 


0.55 


0.0098 


0.0040 


+0.9 






VLBA observed 


on May 09. 


2005 














1.667L 




85.7 


± 


4.4 


0.58 


0.0104 


0.0039 


-1.4 






VLA-B observed 


on Dec 15 


, 2003 














44.94 




334.1 


± 


10.1 


0.56 


0.9 


0.4 


-48.1 






8.460 




359.8 


± 


10.8 


0.36 


1.8 


0.7 


-50.3 






4.860 




304.7 


± 


9.7 


1.89 


3.0 


1.2 


-49.9 






1.425 




176.6 


± 


6.0 


1.53 


9.3 


4.2 


-50.6 







Col. (1) observing frequency. L and R denote left and right circular polarization, respectively; Col. (2) component name. CO and CI denote 
the bright and weaker component, respectively, which were spatially resolved by visibility-based model fitting (Section 3.1); Col. (3) flux 
density; Col. (4) RMS of image noise; Col. (5) FWHM of synthesized beam at major axis; Col. (6) at minor axis; Col. (7) position angle 
of synthesized beam; Col. (8) separation between C0-C1; Col. (9) position angle of CI from CO. 



from 111.5 mJy at 1.4 GHz by the FIRST survey (Becker 
et al. 1995) and 234 mJy at 8.4 GHz by the Jodrell Bank 
VLA Astrometric Survey (Browne et al. 1998) showed an 
inverted spectrum. Simultaneous observations at 2.7, 5, 
and 10.7 GHz (Reich et al. 2000) showed a spectral index, 
a = —0.24 + 0.08. Our VLA flux densities were consistent 
with or less than those of the VLBA; this suggests that 
there is no missing flux in the VLBA images of this object. 

The VLBA flux density at 1.7 GHz decreased to 50% 
in 19 months. The flux density at 15 GHz also decreased 
to 69% in 46 days between the VLBA and VLA obser- 
vations. We believe that these were intrinsic decreases of 
flux density rather than losses of high-frequency coherence 
on the VLA data. The VLA observations were made in 
good weather conditions. The coherence loss should be 
inhibited to the same level as that of the VLBA, since a 
sufficiently short solution interval, 10 s, was adopted for 
phase self-calibration to both the VLBA and VLA data 
(Section 2). 

4. Discussion 

4-1. Very high brightness temperatures 

We can derive brightness temperatures from both the 
VLBA image and the flux variability. Using imaging anal- 
ysis, brightness temperature in the rest frame is calculated 
using (cf. Lahteenmaki et al. 1999) 

1-1-7 r 2 qobs 

T Brcst (img) = I 2 ,r bS /9 , 2 ■ (!) 

where c is the speed of light, fc B is the Boltzmann con- 
stant, S u is the flux density at frequency v 1 9 is the 



angular diameter of a source, and suffixes of "obs" and 
"rest" denote quantities in the observer frame and the 
rest frame, respectively The Dopplcr factor is defined as 
5 = \J\ — (3 2 1 (1 — (3co&4>), where j3 = v/c (v is the source 
speed), and <f> is the "viewing angle" between the direc- 
tion of the source velocity and our line of sight. Given 
an apparent diameter of <0.10 mas for CO in the 15-GHz 
VLBA data (Section 3.1), T Brcst (img) > 5.5 x lO n /5 K. 

Alternatively, with flux variation, a brightness temper- 
ature in the rest frame can be calculated independently 
using (cf. Lahteenmaki et al. 1999) 

T Brest var) > — - — — ° ba , 2) 

5 6 27Tfc B ^obs( A< °bs) 

where D co is the comoving distance, AS 1 "^ is a change of 
observed flux density in a period, of A< b s - In this equa- 
tion, we assume exponential decay as seen in radio out- 
bursts in blazars (Valtaoja et al. 1999). We found that the 
flux density decreases by 31% in 46 days at 15 GHz, and 
by 50% in 19 months at 1.7 GHz (Section 3.3), which im- 
plies that T Brost (var) > 3.3 x 10 13 /c5 3 K and T Brcst (var) > 
1.1 x 10 13 /<5 3 K, respectively. A very high brightness tem- 
perature in this source has also been suggested by Zhou 
et al. (2003) from 1.4-GHz flux variation in ^years time 
scale. These very high brightness temperatures form clear 
evidence for a nonthermal process in the central engine of 
this NLS1. 

4-. 2. Highly relativistic jets 

Based on near equipartition of energies between radi- 
ating particles and a magnetic field, the brightness tem- 
perature of nonthermal radio emission should be within a 
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physically realistic upper limit of ~ 2 x 10 11 K in the rest 
frame (Readhead 1994). Therefore, the Doppler factor 
should be 5 > 2.7-5.5 for the observed brightness temper- 
atures (Section 4.1). The observed radio emissions are 
presumably from highly rclativistic jets and are enhanced 
by Doppler boosting. The jets should approach us with a 
viewing angle of 9 <22° and a speed of (3 > 0.76 for 5 > 2.7 
(Section 4.1). CI, spatially resolved in the VLBA data 
(Section 3.1), may be one of the expanding jet compo- 
nents. 

Because of the very high brightness temperature, syn- 
chrotron self-absorption could be a possible explanation 
for the inverted spectrum as observed. We show an exam- 
ple of spectral model fitting to the observed VLBA total 
fluxes using power-law spectra with a low-frequency cut- 
off due to synchrotron self-absorption: S v = S Vo u 2 ' 5 [l — 
cxp (— Ti/ a °~ 2 - 5 )], where S vo is a scaling constant, r is the 
optical depth, and ao is the intrinsic (not absorbed) spec- 
tral index. Spectral model fitting was performed on the 
assumptions that the observed spectrum of CO is caused 
by a superposition of self-absorbed synchrotron spectra 
of two components, and that the observed spectrum of 
CI can be represented by a single power-law spectrum. 
As a result, CO consists of two convex spectra with peak 
frequencies at -16 GHz ("COA") and -3 GHz ("COB"), 
as shown in Fig. 2. The Doppler boosting can be re- 
sponsible for such high frequencies of the spectral peaks, 
since v b s = [5/{l + z)]is rcst . The intrinsic peak frequen- 
cies could be at several hundreds of MHz. Such high- 
frequency peaked spectra are reminiscent of those of the 
Seyfert galaxy III Zw 2 with two compact radio compo- 
nents and an inverted high-frequency spectrum in total 
flux (Falcke et al. 1999). These components in III Zw 2 
are also highly rclativistic jets (Brunthalcr ct al. 2000). 

4.3. Implications of SDSS J094857. 3+002225 

Our observations showed a slightly inverted spectrum at 
-5 GHz and Rk 2530, which is derived from m B (SDSS) « 
18.89 mag (Zhou et al. 2003) and the VLBA total flux 
density at 5 GHz. Given 5 > 5.5 (Section 4.1), the peak 
frequency of COA, ~16 GHz, should be <2.9 GHz before 
Doppler beaming. In the case where a = —0.6, a flux 
density should be amplified by a factor of >462 from an 
intrinsic spectrum at the same frequency. Therefore, there 
is a possibility that SDSS J094857.3+002225 is a radio- 
quiet object with a steep spectrum in the rest frame, like 
typical NLSls (Moran 2000; Zhou & Wang 2002; Ulvestad 
ct al. 1995). 

Flat or inverted radio spectra have also been found in 
the other radio-loud NLSls, e.g., RGB 0044+193 (Siebert 
et al. 1999), 2E 1640+5345, RXS J16290+4007, and 
RX J16446+2619 (Zhou & Wang 2002). Doppler boost- 
ing may be responsible for their spectra and radio loud- 
ness as well as those of SDSS J094857.3+002225. The ob- 
served optical and X-ray properties and black-hole masses 
of radio-loud NLSls are not significantly different from 
those of radio-quiet ones (e.g., Siebert et al. 1999; Wang 
et al. 2001). From our observations, we suggest that the 
Doppler factor is a crucial parameter determining the ra- 
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Fig. 2. Radio continuum spectra of 

SDSS J094857.3+002225. Symbols representing total 
flux densities are identified in top corner. Filled circles 
particularly describe the flux densities for component CI. 
The solid curve represents the spectral model that was 
determined by spectral model-fitting to the VLBA total 
fluxes using two self-absorbed synchrotron spectra ("COA" 
and "COB" as dashed curves) with an intrinsic spectral index 
of ao = —0.6 and one simple power-law spectrum (CI as a 
dashed line). The spatially single component CO was divided 
into two spectrally separate components, COA and COB 
with peak frequencies of ~16 and 3 GHz, respectively. See 
Section 4.2 for details. 

dio loudness of NLSls: at least some of NLSls could pos- 
sess the ability to generate highly relativistic jets. 

Analogies have been drawn between X-ray binaries and 
AGNs (e.g., Maccarone et al. 2003). NLSls are thought 
to be AGNs with high accretion rates. Although X- 
ray binaries in the high/soft state with accretion rates 
of L/L^dd ~ 0.03-0.1 have no relevance to jets, strong 
radio flares have been seen during the very high state, 
which is thought to be a disk in L/L^dd > 0.1, such as 
GRS 1915+105, with highly relativistic jets (e.g., Fender 
2003). SDSS J094857.3+002225 has the potential to be 
compatible with the very high state, rather than the 
high/soft state. Radio-loud NLSls will be appropriate tar- 
gets to investigate for the very high state on super-massive 
black hole systems. 

5. Summary 

Our VLBA and VLA observations have revealed that 
SDSS J094857.3+002225, which shows the highest known 
radio loudness for an NLS1, has (1) significant variabil- 
ity on time scales of 46 days and 19 months, (2) an in- 
verted spectrum, and (3) multiple components with a 
size <0.10 mas. The derived brightness temperatures 
were significantly higher than the intrinsic limit, implying 
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Doppler-beaming with a Dopplcr factor S > 2.7-5.5 and, 
therefore, relativistic jets. At least some NLSls could be 
able to generate highly relativistic jets. We suggest that 
the Doppler factor is a crucial parameter determining the 
radio loudness of NLSls. On the analogy of X-ray bina- 
ries, SDSS J094857.3+002225 is probably equivalent to 
an AGN in the very high state rather than the high/soft 
state, because of the highly relativistic jets. 

The VLBA and VLA are operated by the National 
Radio Astronomy Observatory (NRAO), a facility of the 
National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. This research 
has made use of NASA Astrophysics Data System (ADS) 
bibliographic services. We thank the referee for sugges- 
tions that helped considerably in clarifying the paper. 
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